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*  Description  of  shear  coaxial  jets 

*  Applications  to  rocket  engines 

*  Survey  of  relevant  previous  works 
—  Experiments 

—  Numerical  simulations 

*  AFRL’s  published  research  in  the  area 

*  Current  challenges 

*  Conclusions 
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Anatomy  of  shear  coaxial  jets 
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From  J.  I.  Rodriguez.  Acoustic  Excitation  of 
Liquid  Fuel  Droplets  and  Coaxial  Jets.  Ph.D. 
dissertation,  UCLA,  2009. 


Dahm  et  al,  JFM,  Vol  241, 1992 
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Important  flow  and  geometry 
parameters 


1 .  AR:  Outer  to  inner  jet  Area  Ratio 

2.  J:  Outer  to  inner  jet  Momentum 
Flux  Ratio,  p0U20/  pJJ2^ 

3.  R:  outer  to  inner  jet  velocity  ratio, 

UJ  Ux 

4.  Re:  Reynolds  number  for  inner  and 
outer  jets,  p]\J] 

5.  S:  outer  to  inner  jet  density  ratio, , 

Pj  Pi 

6.  t/Dp  Inner  jet  lip  thickness/inner 
diameter 

7.  We:  Weber  number  if  surface 
tension  is  present, 

8.  Exit  boundary  layer  profile  -  3 
relevant  boundary  layers 

9.  Thermodynamic  states  (2  phase,  1 
phase) 


1 0.  Transverse  Acoustic  mode  from 
chamber/siren,  f=f(c,  geometry 

St=fDij/Uij 

1 1 .  Acoustic  modes  for  outer  and  inner 
jets 

f~c/2L  -  2  speeds  of  sound  up  and 
downstream 

12.  Wake  resulting  from  inner  post 
thickness 

St=fw/Uch 

13.  Shear  layer  instabilities 
Ste=f0/Uch 

14.  Jet  preferred  modes 

St=fDjj/Ujj 
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Some  applications  in  rockets 


Vulcain  2  from  EADS 
ASTRIUM 

http://cs.astrium.eads.net/sp/lau 

ncher-propulsion/rocket- 

engines/vulcain-2-rocket- 

engine.html 


B.  K.  Wood,  Propulsion  for  the  21st 
Century— RS-68,  AIAA  2002-4324 

Space  Shuttle  Main 
Engine  Injector  face 

Rocket  Propulsion  Elements  (7th  Edition), 

Sutton,  G.  P.;  Biblarz,  O.  2001 


HM  7  from  EADS 
ASTRIUM 

http://cs.astrium.eads.net/sp/lau 

ncher-propulsion/rocket- 

engines/hm7b-rocket- 

engine.html 
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SSME  view  from  nozzle  exit 


6 


shear 


coaxial 


injectors 


Injector 
Face  with 


SSME 

http://www.nasa.gov/images/content/68 

098main_nasm_ent_ssme_hires.jpg 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


What  do  we  need  to  know  about 
shear  coaxial  injectors? 
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The  mixing  process  between  the  inner  and  outer  jets  for  different 
conditions: 

If  applicable,  need  to  characterize  droplet  formation,  droplet 
diameter  distribution  and  vaporization 

Modern  rockets  operate  at  supercritical  pressures  with  respect 
to  the  propellants  therefore  also  need  to  understand  mixing  and 
combustion  beyond  liquid  and  gas  states 

An  injector’s  own  instability  modes 

The  interactions  of  transverse  acoustics  with  injector’s  own 
instability  modes  and  mixing  processes 

Need  to  understand  differences  in  response  to  pressure  and 
velocity  nodes 

What  non-dimensional  numbers  are  most  critical  to  define  mixing 
and  instability  modes 

How  injector  geometry  affects  mixing  and  susceptibility  to  external 
disturbances 
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Consequence  of  incomplete  knowledge  - 

Combustion  Instabilities 
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Injector  head 
molten  after  a 
combustion 
instability 
event 


Courtesy:  U.S.  Rocket  and  Space  Center,  Huntsville,  AL 


Combustion  instability  is  an  unsustainable  cimwth  of 
pressure  anc/  heat  transfer  fluctuations  in  a  rocket  engine 


Irreparable  damage  can  occur  in  <1s 
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Earlier  contributors  to  the  field  1/2 
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Michalke,  1964 

—  Linear  stability  theory  for  inviscid  instability  of  a 
hyperbolic  tangent  velocity  profile  for  free  boundary 
layers 

•  U(y)=0.5[1  +  tanh(y)] 

Chigierand  Beer,  1964 

—  pioneering  experimental  work  on  uniform  density  air 
coaxial  jets 

—  detailed  velocity  and  static  pressure  measurements  of 
fully  developed  turbulent  flow  with  t/D1~1 

—  After  the  two  jets  combine  the  behave  as  an  equivalent 
single  jet 

—  As  R  increases  the  inner  potential  core  decreases 

—  At  a  short  distance  from  jet  exit  each  jet  can  be  treated  as 
separate  jets 


From  Chigier  NA.  and  Beer  JM,  The  Flow  Region  Near  the 
Nozzle  in  Double  Concentric  Jets,  J  of  Basic  Eng,  1964 


•  Crow  and  Champagne,  1971 

—  Single  jet  preferred  mode  “large-scale 
vortex  puffs”,  Std=fd/U~0.3 

—  For  Re:[102-103]  instabilities  go  from 
sinusoid  to  helix  to  axisymmetric  waves 

—  Studied  response  to  inline  acoustics 


From  Crow  SC,  Champagne  FH,  Orderly  structure  in  jet  turbulence,  JFM,  1971,  (48)  547-591 
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Earlier  contributors  to  the  field  2/2 
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Ko  et  al,1 976-1 989 

—  Some  of  the  most  extensive  earlier  work  by  a  single 
group  on  coaxial  jet  characterization 

—  Very  detailed  description  of  near  field  mixing  for  coaxial 
jets.  Many  later  authors  cite  his  nomenclature 

—  Similarity  solutions,  variation  of  St  with  axial  distance 

—  Development  of  vortex  structures  in  the  inner  and  outer 
shear  layers 
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Kwan  and 
Ko,  J. 
Sound  and 
Vibration, 
48  (2),  1976 


1 .  A.  S.  H.  KWAN  and  N.  W.  M.  Ko  1976  Journal  of  Sound  and  Vibration  48,  203-219.  Coherent  structures  in  subsonic  coaxial  jets. 

2.  Kwan  ASH  and  Ko  WM,  J.  Sound  and  Vibration,  48  (2),  1976,  Coherent  structures  in  subsonic  coaxiial  jets. 

3.  N.  W.  M.  KO  and  A.  S.  H.  KWAN  1976  Journal  of  Fluid  Mechanics  73,305332.  The  initial  region  of  subsonic  coaxial  jets. 

4.  Kwan,  A.  S.  H.  and  Ko,  N.  W.  M.,  "The  Initial  Region  of  Subsonic  Coaxial  Jets.  Part  II,"  Journal  of  Fluid  Mechanics,  Vol.  82,  No.  2,  1977,  pp,  273- 
287 

5.  Au,  H.  and  Ko,  N.  W.  M.,  "Coaxial  Jets  of  Different  Mean  Velocity  Ratios.  Part  2,"  Journal  of  Sound  and  Vibration,  Vol.  116,  No.  3,  1987,  pp.  427- 
443. 

H  Au,  “Coaxial  Jets  of  Different  Mean  Velocity  Ratios”,  J  of  Sound  and  Vibration,  Vol  100  No.  2,  1985,  211-232 
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Highlights  from  previous  work  -  70s  to  early 
80s 


Boldman  et  al,  1975 


—  Experimental  and  theoretical  analysis  for  mixing  of  two  air  streams  with  different 
velocities  -  points  out  different  vortex  interactions,  St,~0.2  (Uave) 


-j  h 
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Gutmark  and  Ho,  1983 


Early  analysis  of  what  happens  a 
recirculation  zone  when  two 
streams  separated  by  a  blunt  edge 
are  of  unequal  velocities 


Even  those  the  visualizations  were 
not  very  clear  the  analysis  captures 
the  change  in  structures  that  later 
experiments  will  verify 


Wi 
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From  Boldman,  D.R.,  Brinidh,  P.F.,  and  Goldstein,  M.E. 
"Vortex  Shedding  from  a  Blunt  Trailing  Edge  with  Equal 
and  Unequal  External  Mean  Velocities,"  JFM,  (75)  4, 
1976,  721-735 


—  Collects  previous  results  on  jet  preferred  mode,  Std  has  a  range  from  ~0.24-0.64 
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Highlights  from  previous  work  -  early  90s 
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Dahm  et  al,  1992 

—  Studied  different  instabilities  in  the  inner  and  outer  shear  layers  plus  effect  of  absolute 
velocity  and  R 


Dahm  et  al,  JFM  (241)  1992 


R  increases  > 

Wicker  and  Eaton,  1994 

—  Forces  air  inner  and  outer  jets  independently  -  observes  vortex  growth 

Wicker  and  Eaton.  AIAA  J.  (32)  No.3.1994 
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Highlights  from  previous  work  90s  -2000s 


13 


Buresti,  Talamelli,  Petagna  (1994, 1998) 

—  Air  jets,  same  density,  top-hat  profile, 
Stdo=fdo/Uoj~0.3  to  1,  function  of  x/Dj 


>T|pir*  13  VamilJO*  Svwtol  rainiwr  Si,  nil  ijtmjrr,n|. 
jiiKi  wtodij. 


Buresti  et  al,  Exp  Thermal  and  Fluid  Science, 
1994,  9:135-146 


Balarac,  da  Silva,  Metais  et  al  (2003,  2007) 

—  DNS  analysis  of  coaxial  jets  -  same  density,  top-hat  profiles 

—  Consider  two  shear  layers,  study  effect  of  R 

—  Consider  axisymmetric  and  azimuthal  excitation 


I  ipurt  4  hhLir-.Lim  ui  LUCd-ixi  uf  a  nLiJz  fr*  :  n  Tir  DNSIO  uk.  •  lI  dl  iH  .UiUaJ  JlLaf*:  iltl  ft  ■  m  Us 
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rri.  -T*  —  -  *-»  |pn  =1-  imp-ML 


Segalini,  Talamelli,  et  al  (2006,  2011) 

—  Air  jets,  same  density,  top-hat  profile,  Stb(|ip)=fb(lip)/Uaverage 


Segalini  et 
al,  Phys  of 
Fluids  (23), 
2011 
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Highlights  from  previous  work,  late 
2000s 


*  Birbaud,  Ducruix,  Durox,  Candel  (2006-2007) 

—  Single  air  jets,  low  Re,  laminar,  top-hat  profile,  subjected  to  acoustic  modulation 
—  Systematic  study  of  effect  of  modulation  in  terms  of  Std,  Ste 
—  Need  similar  study  for  shear  coaxial  jets 


Highlights  from  previous  work  late 
2000s  to  today 


*  Richecoeur,  Scouflaire,  Ducruix,  Candel  (2006) 

—  Forced  transverse  acoustic  excitation  of  flames 


*  Kastengren  etal2013 

—  X-ray  radiography  used  to  quantitatively  examine  shear  coaxial  jet  injectors  First 

quantitative  measurement  of  the  mass  density  in  the  near  injector  region  of  this  injector 
type.  Centerline  EPL  profiles  showed  same  general  trend  with  J  reported  in  studies  of  the 
dark  core  length.  Shortening  of  the  core  region  with  increasing  J. 
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AFRL’s  previously  published  research  on  shear 
coaxial  jets:  non-reacting  studies 
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•  Study  the  mixing  process  and  response  to  external  acoustic 
perturbations  for  shear-coaxial  injectors  typically  used  in 
cryogenic  engines 

•  Investigate  influence  of  injector  geometry 

•  Vary  the  outer-to-inner  jet  momentum  flux  ratio,  Jfor  P<  1  and 
Pr>  1 


3  A  MPa  (493  psi) 


critical  ,^2 


p  - ^nu^nu^ 

'r  p 

critical, 


•  Characterize  mixing  using  dark-core  length  measurements 

•  Apply  proper  orthogonal  decomposition  of  high-speed  image  pixel 
intensity  fluctuations  to  extract  spatial  and  temporal  characteristics 
of  flow  structures 
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Image  of  experimental  facility 
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From  Teshome,  S.,  Droplet 
Combustion  and  Non-Reactive 
Shear-Coaxial  Jets  with  Transverse 
Acoustic  Excitation,  Ph.D.  thesis, 
2012,  UCLA  (Teshome  20 12 A) 


Schematic  of  experimental  facility 


18 


_ 


/ 

Piezo-Siren 


GN2 

(Ambient  T) 


LM 


Out 


ner 


Outer 

Jet 


0 


Xenon 

^rc-Lamp 


Main 

Pressure 

Chamber 


From  Teshome,  S.,  2012A 
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Density  gradients  around  critical 
point  of  nitrogen 
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Injector 

Di 

D2 

tID \ 

AJA, 

IID \ 

LAR-thick 

0.51 
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2.42 

3.18 

1.05 

12.9 

0.5 

SAR-thin 

1.40 

1.65 

2.44 

3.94 

0.09 

1.6 

0,0.5 

SAR-thick 

1.47 

3.96 

4.70 

6.35 

0.84 

2.9 

0.5 

LAR-thin 

0.70 

0.89 

2.44 

3.94 

0.13 

10.6 

-0.2 

SAR,  LAR  Small,  Large  Area 
Ratio 

Thick,  Thin  Post  thickness 

From  S.  Teshome,  I  .A.  Leyva,  J.  Rodriguez,  D. 

Talley,  Geometry  Effects  on  Steady  and 
Acoustically  Forced  Shear-Coaxial  Jet  Sprays, 

ICLASS,  2012 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Dark-core  length  measurement 
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•  First  raw  grayscale  images  were  converted  to  binary  images 

•  A  contour  was  drawn  around  the  “dark-column”  in  the  binary  image 

•  Axial  length  of  the  dark-column  measured  and  defined  as  the  Dark-Core  Length,  L 

Grayscale  (8-bit)  Binary  Image  Raw  Image  w / 


From  Teshome  2012A 
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Break-up  regimes  for  shear  coaxial 
jets 
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We 


From  Lasheras,  J.  C.,  and  Hopfinger,  E.  J.,  “Liquid  Jet  Instability  and  Atomization  in  a 
Coaxial  Gas  Stream,”  Annual  Rev.  Fluid  Mech.,  Vol.  32,  2000,  pp.  275-308. 
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Different  correlations  for  inner  jet 
break-up  length 


23 


From  D.W.  Davis.  On  the 
Behavior  of  a  Shear-Coaxial  Jet, 
Spanning  Sub-  to  Supercritical 
Pressures,  with  and  without  an 
Externally  Imposed  Transverse 
Acoustic  Field.  Ph.d. 
dissertation,  Pennsylvania  State 
University,  2006. 
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State  of  the  art  for  mixing  in  shear 
coaxial  jets  -  2007 
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From:  I.  Leyva,  B.  Chehroudi,  D.  Talley,  Dark  Core  Analysis  of  Coaxial  Injectors  at  Sub-,  Near-,  and  Supercritical  Conditions 
in  a  Transverse  Acoustic  Field,  AIAA  2007-5456 
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Differences  between  2-phase  and  1  - 
phase  coaxial  jets 
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From  J.  I.  Rodriguez.  Acoustic  Excitation  of  Liquid  Fuel  Droplets  and  Coaxial  Jets.  Ph.D.  dissertation,  UCLA,  2009. 
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Effect  of  acoustics  on  a  single 
injector  geometry 
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Baseline  dark-core  lengths  for  4 
geometries  for  Pr  =  0.44 
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□  SAR-thin  (Rodriguez) 

+  SAR-thin  (Graham  et  al.) 
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Fastest  Mixing  for  SAR-thick  -  Slowest  Mixing  for  LAR-thin 
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Baseline  normalized  DCL,  L  JD1  - 
Pr  =  0.44 


From  Teshome  2012A 
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Baseline  flows  at  Pr  =  0.44 


LAR-thin  Injector 


SAR-thick  Injector 


Baseline  normalized  DCL,  L  JD1  - 
Pr-  1.05 


LAR-thin 

o 

LAR-thick 

® 

SAR-thin 
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Baseline  flows  at  Pr  =  1 .05 


LAR-thin  Injector 


J=  0.1  0.5 


SAR-thick  Injector 


Baseline  normalized  DCL,  LBID1 
collapse 


Has  characterized  mixing  for  shear  coaxial  injectors  taking  into  account 

geometry  and  flow  conditions! 
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ACOUSTIC  EXCITATION 

ANALYSIS 
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Acoustic  field  set-up:  pressure 
antinode 


*  Pressure  antinode  (PAN)  -  condition  of  maximum  pressure  perturbation 
in  the  acoustic  field 

*  Piezo-sirens  forced  in-phase 

*  Superposition  of  quasi-1  D  acoustic  waves  traveling  in  opposite 
directions  =>  PAN  at  the  jet  location  (geometric  center  of  test  section) 


Standing  Acoustic  Pressure 
Waves 


'Piezo-Siren 


Time  Trace  of 
Pressure 
Fluctuations 
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Proper  orthogonal  decomposition 
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•  Proper  Orthogonal  Decomposition  (POD  was  used  for  extracting  dominant 
dynamical  processes  embedded  in  high-speed  images. 

•  A  time-resolved  set  of  images  A(x,t )  can  be  represented  as  a  linear  combination 
of  orthonormal  basis  functions  (f)k  (aka  proper  orthogonal  modes)1’2’3  : 


M 

A(x,  /)  =  £>*  (t)A  (x) 

k= 1 

where  ak(t)  are  time  dependent  orthonormal  amplitude  coefficients  and  M  is  the 
number  of  modes 


•  Main  idea:  POD  modal  amplitudes  capture  the  maximum  possible  “energy”  in 
an  average  sense4,  i.e., 

X(a,(0«,(0}>Z(W*(0) 

k  k 

where  bk(t)  are  the  temporal  coefficients  of  a  decomposition  with  respect  to  an 
arbitrary  orthonormal  basis  y/k. 
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2  Arienti,  M,  and  Soteriou,  M.C..  Phys.  Fluids  21,  112104  (2009) 
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Construction  of  data  set 
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•  First,  form  a  row  vector  consisting  of  all  pixel  intensity  values  of  each 
snapshot  image  (with  resolution  of  n  rows  by  m  columns)  in  order  of  increasing 
columns,  then  increasing  rows 


Image 


columns 


t 

n 

rows 

i 


•  Then,  combine  all  such  row 
vectors  for  N  sequences  of  image 
frames  resulting  in  a  matrix  A 
consisting  of  N  rows  by  P  =  n*  m 
columns  of  intensity  values. 


A  = 


\7 

▲ 

_ P  =  n  x  m 

pixel  intensities 

N  time 
steps 


*■ 
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Orthogonal  decomposition 
technique 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 

•  SVD  preferred  since 

1 .  Applicable  to  non-square  matrices  (most  likely  the  case) 

2.  Decomposition  matrices  are  orthogonal 

3.  Subroutine  readily  available  in  MATLAB® 

•  Subtracted  temporal  mean  of  A  =>  matrix  of  intensity  fluctuations  A 


•  Applied  SVD 


=>  A  =  USVT  =  QVt 

Columns  of  Q  ~  ak(t)  contain  temporal 
information 


N  x  N 


Orthogonal  Matrix  of  Left 
Singular  Column  Vectors  of 
A 


Columns  of  V  ~  <f>k(x <)  contain  spatial 

*>  N  x  p  information 

Orthogonal  Matrix  of  Right  Singular 

Column  Vectors  of  A  proper 
orthogonal  modes  (POM) 


A/x  N 

Diagonal  Matrix 
of 

Singular  Values  From  Teshome  2012A 
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Results  -  LAR-thin,  Pr  =  0.44,  J  =  5.2 


LAR-thin 


Snapshot  Average 


POM  1  POM  2 


Amplitude  information 
contained  in  singular  values 


Antisymmetric 
Structures  Identified  with 
Characteristic  Frequencies 


10  10 
Mode 


Power  Spectral 
Densities  (PSD) 
of  Temporal 
Coefficients  of 
POMs  1  and  2 
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Cross-Power  Spectral  Density 
(CPSD) 

- - - -  LAR-thin 


CPSD  yields  the  FFT  of  the  cross-correlation  of  the  temporal 

coefficients*  ^4cov(a  ,a  )  .  , 

CPSD  =  £ -  "  -1  eiak 


k= 0 


<y  a 

a v  a , 


•  Magnitude  and  phase  plots  used  to  determine  existence  of  propagating 
LAR.  Pr  =  0.44.  J=  5.2  structures  pAN 


*Arienti,  M,  and 
Soteriou,  M.C.. 

Phys.  Fluids  21, 
112104  (2009) 


0  3  6  9 

/(kHz)  From  Teshome  2012A 
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Results  -  LAR-thin  Pr=0.44 


Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities  for 
all  J 


J—  2.1 


/=  20 


Baseline 


PAN 


Baseline 


PAN 


/(kHz) 


106 


/(kHz) 


L 


Geometry  with  least  efficient  mixing  -  less  responsive  to  acoustics 

at  high  J 
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CPSD  Magnitude 


Results  -  SAR-thin  Pr=0.44 


Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 


J=  2.0 


J=u 


Baseline 


PAN 


Baseline 


PAN 


(fF  =  2.97  kHz) 


(fF  =  2.90  kHz) 


Results  -  SAR-thick,  Pr=0.44 


SAR-thick 


Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 


J=  2.1 


j=  21 


Baseline 


PAN 

(fF  =  3.07  kHz) 


Baseline 


PAN 

(fF  =  3.11  kHz) 
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Geometry  with  most  efficient  mixing  -  more  responsive  to 

acoustics  at  high  J 
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Collaboration  with  Caltech 
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*  Simplification  of  axisymmetric  coaxial  jets  to 

VR=i  a  2D  setup 

•  Investigated  effect  of  velocity  ratio  on  St 
number  on  a  water  tunnel  with  PIV 


VR=0.55 


VR=0.33 


Velocity  ratio  vs  Strouhal  number 


Differences  in  recirculation  zone  evident  as  VR 
decreases  -  St  decreases  with  VR 


From:  Split  Stream  Flow  Past  a  Blunt  Trailing  Edge  with  Application  to 
Combustion  Instabilities,  V.  Tian,  B.  McKeon,  I.  Leyva,  AIAA  20 12-3807 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Collaboration  with  Ecole  Centrale  Paris:  comparison 
between  CFD  and  experiments  for  shear  coaxial  jets 


CENTRALE 

*  A  k  I  B 
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Velocity  Fields  for  Different  J  values 
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Longitudinal  slice  of  axial  velocity.  White:  minimum;  black:  maximum. 

Dark  line  indicates  iso-contour  of  zero  axial  velocity.  Top:  J  =  1.1;  middle:  J 

=  3.0;  bottom:  J  =  9.3. 

From  Experiments  and  numerical  simulation  of  mixing  under  supercritical  conditions,  T. 
Schmitt,  J.  Rodriguez,  I.  A.  Leyva,  and  S.  Candeh  Physics  of  Fluids  24,  2012 
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Top:  centerline  profile  of  density;  bottom: 
centerline  profile  of  velocity.  Dark  line:  J  = 
1.1;  dash  point  line:  J  =  3.0;  dashed  line:  J  = 
9.3. 


Fundamental  frequencies  for 
baseline  conditions 


Case  “N6” 
(J=3.05) 
Density 
distribution 
(white:  60 
kg/m3; 
black:  410 
kg  /m3; 
logarithmic 
scale). 
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Case  “N8” 

( J  =  9.3) 
Density 
distribution 
(white:  60  kg/ 
m3;  black: 
maximum ; 
logarithmic 
scale). 
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Found  relevant  St  numbers  for  our  configurations 
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•  Research  on  shear  coaxial  jets  goes  back  to  60s 

•  First  researchers  focused  on  the  development  of  the  jets  for  different 
variables:  R,  geometry 

*  Defined  the  different  mixing  regions  in  a  shear  coaxial  injector  -  inner  and  outer  potential 
cores 

*  Made  a  lot  of  very  rich  visualizations  and  measurements 

*  Found  functional  relations  between  R  and  the  potential  core 

9  A  lot  of  previous  work  investigates  the  relevant  St  for  the  inner  and  outer 
jets  and  what  happens  when  jets  are  modulated  close  to  their  preferred 
modes 

9  In  the  last  ten  years  a  lot  of  progress  has  been  made  on  LES  and  DNS  to 
better  understand  the  mixing  processes  in  shear  coaxial  jets 

9  AFRL  has  published  results  on  four  geometries  with  the  primary  goal  of 
understanding  mixing  in  unforced  and  acoustically  driven  shear  coaxial 
jets 

*  The  mixing  process  is  measured  qualitatively  by  the  dark  core  which  is  related  to  the 
potential  core  for  both  unforced  and  driven  conditions 
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•  J  is  not  the  only  parameter  to  characterize  mixing  in  shear  coaxial  jets  - 
area  ratio,  inner  jet  post  thickness  and  density  also  matter 

•  Found  functional  correlations  for  dark  core  lengths 

9  Proper  orthogonal  decomposition  of  high-speed  image  intensity 

fluctuation  data  revealed  key  spatial  and  temporal  characteristics  of  flow 
structures 

9  Based  on  four  geometries,  mixing  is  enhanced  with  increasing  J  and  t/DI 
but  less  efficient  with  increasing  Ao/Ai 

•  Helical  instabilities  were  observed  for  Pr=0.44  for  LAR-thin,  SAR-thick, 
and  SAR-thin  geometries  at  large  J  for  baseline  conditions 

•  Geometry  plays  a  role  for  the  response  of  the  jet  to  acoustics,  at  Pr=0.44 


PAN: 


9  LAR-thin:  The  flow  became  less  sensitive  to  applied  PAN  forcing  with 
increasing  J 

•  SAR-thin:  Flow  responds  to  forcing  acoustic  frequency  at  all  J 

9  SAR-thick:  Recirculation  zone  dampens  the  effect  of  PAN  forcing  at 
lower  J,  but  at  higher  J  (>10)  measurable  effect 
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Baseline  flows:  LAR-thin  injector 
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Pr  =  0.44 


J=0.1 


Pr  =  1 .05 


□  □ 


Baseline  flows:  SAR-thick  injector 
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Pr  -  0.44 


Image  interpretation  key 


Baseline: 

Acoustics 

off 


0° 

Acoustics  on 
PAN,  VN 

“Pressure  coupled” 


180° 

PN,  VAN  " 
Largest  difference 

O 

expected  from  0 
“Velocity  coupled” 


pressure  =  fixed 
J  (p0u02/pjUj2)  =  fixed 

PN  -  pressure  node  -  Min 
PAN  -  pressure  antinode  -  Max 
VN  -  velocity  node 
VAN  -  velocity  antinode 


From  J.  I.  Rodriguez.  Acoustic  Excitation  of  Liquid  Fuel  Droplets  and  Coaxial  Jets.  Ph.D.  dissertation,  UCLA,  2009. 
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Matrix  Decomposition 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 


•  SVD  Subroutine  readily  available  in  MATLAB® 


Prior  to  matrix  decomposition,  the  temporal  mean  of  A  was  subtracted  resulting  in  a  matrix  of  intensity 
fluctuations  A,  i.e.,  _ 

for  /  =  1...N, 


a«=a»-^Za. 


j  =  1...M 


•  Application  of  SVD  on  A  gives  two  orthogonal  matrices  U  (NxN)  and  V  (MxM),  and  a  diagonal  matrix  S 
(NxM)  of  singular  values  in  increasing  order  of  magnitude 


'T'  T' 

A  =  USV  =  QV 

•  Thus,  a  time-resolved  set  of  images  intensity  fluctuations  A(x,f)  can  be  represented  as  a 
linear  combination  of  orthonormal  basis  functions  (/)k  such  that 

M 

A  (x,t)  =  ^ak(t)</>k(x) 

k-\ 


where  ak(t )  are  time  dependent  orthonormal  amplitude  coefficients  and  ^(x)  are  the  proper 
orthogonal  modes  of  A. 


•Equivalence:  columns  of  Q  ~  ak(t),  columns  of  V  ~  ^A(x) 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Collaboration  with  ECP:  grid  and  mesh 


CENTRAL E 

*  tl  1  I  E 


3D  visualization  of  reservoir 


V?-Sifl(2Tf  f  t  + 


Longitudinal  cut  of  domain  with  the  BC's 


Mesh  detail  near  the  injector 


CFD  Injector 


Solver:  AVBP  -  state  of  the  art  LES  code 

*  With  real  fluid  properties  to  tackle  supercritical  fluids 

Mesh:  2  100  000  nodes/10  000  000  tetrahedra 

*  Highly  refined  near  the  injector  (0.032  mm  on  a 
distance  of  10  inner  jet  diameters) 

CPU  hours  on  Europe  Supercomputer 
Center:  100,000 


From  Experiments  and  numerical  simulation  of  mixing  under  supercritical  conditions,  T. 
Schmitt,  J.  Rodriguez,  I.  A.  Leyva,  and  S.  CandeL  Physics  of  Fluids  24,  2012 
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Acoustic  forcing  results 
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T.  Schmitt, 
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Rodriguez, 
I.  A.  Leyva, 
and  S. 
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2012 


Instantaneous  images  of  the  simulation  of  an  acoustic  case  with  the  injector  at  a  pressure  antinode  for  a  J  =  3.0 
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Subcritical  low  J 


55 


Chamber  =  1-5  MPa  (Pr  =  0.44) 


^  J  =  0.089  with  thin  inner  post  injector 


Negligible  effect  vs.  violent  destruction  of  the  jet 
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Subcritical  Moderate  J 
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Chamber  =  1-5  MPa  (Pr  =  0.44) 


J  =  2.6  with  thick  inner  post  injector  J  =  2  0  with  thin  inner  post  injector 

Bending  vs.  vortical  structures  with  fine  atomization 
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Subcritical  large  J 
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Chamber  =  1-5  MPa  (Pr  =  0.44) 
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J  =  23  with  thick  inner  post  injector  J  =  18  with  thin  inner  post  injector 

Two  different  mixing  mechanisms  for  similar  large  J  values 
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Nearcritical  large  J 
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Chamber  =  3-6  MPa  (Pr  =  1.1) 


J  =  9.3  with  thick  inner  post  injector  J  =  9.4  with  thin  inner  post  injector 

Longer  dark  core  lengths  and  hence  visible  effect  of  acoustics 

for  thin  inner  post  geometry 
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Supercritical  moderate  J 
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Chamber  =  5-0  MPa  (Pr  =  1.5) 


J  =  2.4  with  thick  inner  post  injector 


J  =  2.6  with  thin  inner  post  injector 


Large  dark  core  lengths  for  the  thin  post  geometry  prevent  us 
from  observing  acoustic  effects,  if  any 
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Chronological  progression  (only  coaxial 
results  are  summarized  in  what  follows) 


*  Single  jets,  no  coaxial  flow 

—  Davis  et.  al.  (Ph.D.  thesis)  -  single  jets,  no  coaxial  flow 

*  Coaxial  jets 

—  Davis  et  al.  (Ph.D.  thesis)  -  LAR  thick 
—  Leyva  et.al.  -  LAR_thick 

—  Rodriguez  et.  al.  (Ph.D.  thesis)  -  LAR_thick,  SAR_thin 

—  Graham  et.  al.  -  SAR_thin,  two  recesses 

—  Teshome  et.  al.  (Ph.D.  thesis,  expected  April  2012)  - 
complete  all  four  geometries 

*  Also  complete  modal  analysis  of  earlier  geometries 
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•  Mixing  for  shear  coaxial  jets  has  other  major  variables  other  than  momentum 
flux  ratio  (J) 

—  Ratio  of  inner  jet  temperatures  is  another  important  variable 
—  Geometry  -  Area  ratio  and  lip  thickness  also  affect  mixing 

•  For  LARthick 

—  Found  bending  mode  with  largest  effect  at  velocity  antinodes 

—  The  reduction  on  the  dark  core  length  was  greatest  for  a  medium  J  range 

—  For  near  critical  pressures,  the  collaboration  with  ECP  determined  relevant  St  for  our  injector 
configuration  and  was  able  to  capture  qualitative  behavior  of  natural  and  excited  jets 

•  For  SARthin 

—  Did  not  see  bending  mode  for  conditions  studied 

—  Saw  vortex  roll-up  and  puffing  occurring  over  entire  J  range  (0.09-21)  tested  - 
—  PAN  forcing  produced  symmetric  flow  structures  regardless  of  J 


—  Spectral  plots  showed  strong  response  to  PAN  forcing  at  low  and  high  J 

•  For  LARthin 


—  Sees  both  bending  and  vortex  roll-up  modes  depending  on  the  acoustic  frequency 

—  PAN  forcing  at  low  J  produced  symmetric  flow  structures,  while  at  higher  J,  influence  of  forcing 
subsided 

I  magnitude  plots  showed  decreasing  influence  of  PAN  forcing  with  increasing  J 
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